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a.  Principle  Investigator: 

Susan  B.  Sinnott 

The  University  of  Kentucky 

Department  of  Chemical  &  Materials  Engineering 

177  Anderson  Hall 

Lexington,  KY  40506-0046 

b.  Current  telephone  number:  (606)257-5857 

c.  Cognizant  ONR  Scientific  Officer: 

Dr.  P.  P.  Schmidt 
Office  of  Naval  Research 
Ballston  Tower  One 
800  North  Quincy  Street 
Arlington,  Virginia  22217-5660 

d.  Brief  description  of project: 

The  goal  of  the  project  is  to  gain  increased  understanding  of  the  atomic-scale  chemical 
interactions  responsible  for  macroscopic  friction  and  wear  when  two  surfaces  are  brought  into  sliding 
contact  with  each  other.  Such  insight  is  crucial  to  the  development  of  improved  solid-state  lubricants 
and,  consequently,  improved  machines  of  technological  and  military  importance.  The  approach  is 
molecular  dynamics  simulations  to  provide  detailed  information  about  the  tribochemistry  that  occurs 
when  films  of  carbon  solid-state  lubricants,  such  as  graphite  and  diamond-like  carbon,  are  placed 
between  metal  surfaces  in  sliding  contact  with  each  other.  To  facilitate  realistic  interactions  between 
the  carbon  and  metal  atoms,  a  combined  many-body  carbon  -metal  (MBCM)  potential  must  be 
constructed  that  allows  one  to  realistically  characterize  systems  composed  of  thousands  of  carbon 
and  metal  atoms  on  a  standard  workstation. 

e.  Significant  results  during  last  year: 

A  combined  many-body  carbon-metal  (MBCM)  potential  suitable  for  use  in  large-scale 
atomistic  simulations  has  been  constructed.  The  combined  potential  reduces  to  the  reactive 
empirical  bond  order  (REBO)  form  of  Brenner  et  al.  for  carbon-carbon  and  carbon-hydrogen 
interactions,  and  to  the  corrected  effective-medium  (CEM)  form  for  metal-metal  interactions.  The 
metal-carbon  interactions  are  characterized  by  an  expression  of  the  effective-medium  type: 
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E  =  y)[AF(n,)  +  —  y]0(r..)] 

Cohesive  4^  *■  ''  i-  9  —  ‘J 
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(1) 


where  0  is  a  pair  additive  term  that  accounts  for  the  attractive  core-electron  interactions,  and  AF 
is  the  energy  of  "embedding"  each  into  jellium  of  density  n^  (i.e.,  the  energy  difference  between  the 
atom-in-jellium  system  and  the  isolated  atom  and  jellium  systems).  The  jellium  density  is  obtained 
from  the  sum  of  the  overlaps  of  the  atomic  densities. 


The  first  systems  to  be  examined  as  test  cases  are  CH^  (x=0-3)  chemisorbed  on  nickel  surfaces. 
These  systems  were  chosen  because  1)  when  a  metal  surface  is  in  sliding  contact  with  a  carbon 
lubricant,  it  is  expected  that  parts  of  the  lubricant  will  be  transferred  to  the  metal  surface;  2)  there 
is  ample  first  principles  data  available  with  which  to  compare  the  test  results;  and  3)  nickel  substrates 
are  used  to  grow  diamond  films  by  chemical  vapor  deposition  techniques,  a  process  that  could  be 
studied  in  molecular  simulations  with  this  potential.  Some  preliminary  results  are  shown  below: 

(A)  Cohesive  Energy  (eV) 

this  work  first  principles®  this  work  first  principles® 


C-Ni  diatomic 

1.68 

1.63'’ 

-3.49 

-4.93'’ 

C/Ni(100)  [h] 

1.91 

1.85® 

-5.84 

-6.51® 

C/Ni(100)  [a] 

1.73 

- 

-2.55 

CH/Tli(lll)[h] 

1.89 

2.30 

-3.89 

-3.00 

CH/Ni(lll)[b] 

1.83 

2.36 

-2.87 

-2.90 

CH/Ni(lll)  [a] 

1.72 

1.74 

-1.54 

-1.69 

CH,/Ni(lll)  [h] 

1.89 

2.31 

-2.85 

-2.80 

CH2/Ni(lll)  [b] 

1.83 

2.29 

-1.82 

-2.72 

CH2/Ni(lll)  [a] 

1.73 

2.36 

-0.50 

-1.58 

CH3/Ni(lll)[h] 

1.89 

2.59 

-2.76 

-1.68 

CH3/Ni(lll)[b] 

1.83 

2.65 

-1.73 

-1.54 

CH3/Ni(lll)[a] 

1.73 

2.03 

-0.40 

-1.47 

®  From  H.  Yang  and  J.  L.  Whitten,  Surf.  Sci.  255, 193  (1991)  unless  otherwise  indicated. 
”  First  principles  LDA  calculations  from  B.  I.  Dunlap  (private  communication) 

'  From  P.  E.  M.  Siegbahn  and  I.  Panas,  Surf.  Sci.  240,  37  (1990). 
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[h]=hollow  site;  [b]=bridge  site;  [a]=atop  site 

The  table  shows  that  the  MBCM  potential  is  in  good  qualitative  agreement  with  available  first 
principles  data.  Most  of  the  trends  are  captured  correctly,  and  the  potential  could  be  used  to  s^dy 
the  general  tribochemistry  of  a  metallic  material  interacting  with  a  covalently  bonded  material.  It 
is  also  expected  that  the  potential  will  open  up  to  large  scale  computing  many  other  important 
problems  in  materials  and  surface  science,  such  as  corrosion  and  heterogeneous  catalysis.  However, 
some  refining  will  be  necessary  to  obtain  improved  quantitative  agreement  with  the  first  prmciples 
data  and  to  allow  charge  transfer  between  the  atoms. 

f.  Summary  of plans  for  FY  96-97 

Plans  for  FY  96-97  will  be  to  refine  the  MBCM  potential  to  improve  the  agreement  with  the  first 
principles  data  The  first  task  will  be  to  construct  a  new  embedding  energy  fit  for  carbon  (the  first 
term  in  Eq.  (1)).  This  is  currently  fit  to  carbon-carbon  molecular  and  solid-state  data.  A  new  fit  to 
first  principles  data  for  carbon-nickel  interactions  will  lead  to  an  improved  charactenzation  of  the 
carbon-metal  bonds  at  various  electron  density  overlaps.  The  second  task  will  be  to  optimize  the 
bond  angles  for  carbon-carbon  bonds  when  one  of  the  carbon  atoms  is  connected  to  a  metal  atom. 
The  data  to  fit  these  angles  will  also  be  obtained  from  first  principles  calculations. 

Finally,  because  the  chemical  bonding  will  change  drastically  in  the  course  of  the  simulations, 
a  way  of  determining  environment-dependant  partial  atomic  charges  and  charge  transfer  will  be 
necessary  A  formalism  recently  introduced  for  describing  electrostatic  forces  in  alumina  will  be 
adapted  and  incorporated  into  the  MBCM  expression.  This  is  an  approach  that  determines  the  net 
charges  on  each  atom  from  an  environment-dependant  electronegativity  term,  which  is  calculated 
from  the  atomic  electron  densities  associated  with  each  atom.  These  atomic  electronegativites  will 
be  fit  to  reproduce  the  partial  charges  on  the  atoms  predicted  from  first  principles  calculations. 

g.  Graduate  students  and  postdoctoral  associates  associated  with  the  project: 

At  the  moment,  the  Pi’s  group  consists  of  one  PhD  graduate  student,  Mr.  Lifeng  Qi,  and  one 
Masters  graduate  student,  Mr.  Matthew  Flanagan.  It  is  expected  that  Mr.  Qi  will  be  heavily 
associated  with  the  above  project  in  the  course  of  completing  his  PhD  under  the  direction  of  the  PI. 
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Explanatory  material  for  view  graph  information 


Energies  and  Configurations  of  Hydrocarbons  Adsorbed  on  Ni(lOO)  and  Ni(lll) 

A  combined  many-body  carbon-metal  (MBCM)  potential  suitable  for  use  in  large-scale 
atomistic  simulations  has  been  constructed.  The  combined  potential  reduces  to  the  reactive 
empirical  bond  order  (REBO)  form  of  Brenner  et  al.  for  carbon-carbon  and  carbon-hydrogen 
interactions,  and  to  the  corrected  effective-medium  (CEM)  form  for  metal-metal  interactions.  The 
metal-carbon  interactions  are  characterized  by  an  expression  of  the  effective-medium  type.  The 
motivation  behind  the  construction  of  the  potential  is  to  examine  and  characterize  the  atomic-scale 
phenomena  responsible  for  the  macroscopically  observed  friction  and  wear  of  metal  surfaces  in 
sliding  contact  with  each  other  and  the  role  of  carbon-based  solid  lubricants  in  this  process. 

The  view  graph  shows  preliminary  test  results  of  sample  systems  of  CH^  (x=0-3)  chemisorbed 
on  nickel  surfaces.  These  systems  were  chosen  because  1)  when  a  metal  surface  is  in  sliding  contact 
with  a  carbon  lubricant,  it  is  expected  that  parts  of  the  lubricant  will  be  transferred  to  the  metal 
surface;  2)  there  is  ample  first  principles  data  available  with  which  to  compare  the  test  results;  and 
3)  nickel  substrates  are  used  to  grow  diamond  films  by  chemical  vapor  deposition  techniques,  a 
process  that  could  be  studied  in  molecular  simulations  with  this  potential. 

The  first  principles  data  was  obtained  from  the  following  sources: 

•  CH,/Ni(l  1 1)  (x=l-3)  from  Yang  and  Whitten,  Surf.  Sci.  255, 193  (1991) 

•  carbon-nickel  diatomic  from  B.  I.  Dunlap  (private  communication) 

•  C/Ni(100)  from  Siegbhn  and  Panas,  Surf.  Sci.  240,  37  (1990). 
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Results  From  Realistic  Empirical  Carbon-Metal  Potential 
CH,  (x=0-3)  Chemisorbed  on  Ni(lOO)  and  Ni(lll) 
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